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Abstract. Context : The minimum mass for star formation is a critical parameter with profound astrophysical, cosmological 
and anthropic consequences. Aims. Our first aim is to calculate the minimum mass for Primary Fragmentation in a variety of 
potential star-formation scenarios, i.e. (a) hierarchical fragmentation of a 3-D medium; (b) one-shot, 2-D fragmentation of a 
shock-compressed layer; (c) fragmentation of a circumstellar disc. By Primary Fragmentation we mean fragmentation facilitated 
by efficient radiative cooling. Our second aim is to evaluate the role of H, dissociation in facilitating Secondary Fragmentation 
and thereby producing close, low-mass binaries. Methods : We use power-law fits to the constitutive physics, a one-zone 
model for condensing fragments, and the diffusion approximation for radiative transport in the optically thick limit, in order 
to formulate simple analytic e.sXimdiXe.&. Re suits : (i) For contemporary, local star formation, the minimum mass for Primary 
Fragmentation is in the range 0.001 to 0.004 , irrespective of the star-formation scenario considered. This result is remarkable 
since, both the condition for gravitational instability, and the radiation transport regime operating in a minimum-mass fragment, 
are different in the different scenarios. (ii) Circumstellar discs are only able to radiate fast enough to undergo Primary 
Fragmentation in their cool outer parts {R 100 AU). Therefore brown dwarfs should have difficulty forming by Primary 
Fragmentation at S 30 AU, explaining the Brown Dwarf Desert. Conversely, Primary Fragmentation at S 100 AU could 
be the source of brown dwarfs in wide orbits about Sun-like stars, and could explain why massive discs extending beyond this 
radius are rarely seen, (iii) H, dissociation can lead to collapse and Secondary Fragmentation, thereby converting primary 
fragments into close, low-mass binaries, with semi-major axes a ~ 5 AU (^system /O- 1 ), in good agreement with observation; 
in this circumstance, the minimum mass for Primary Fragmentation becomes a minimum system mass, rather than a minimum 
stellar mass, (iv) Any primary fragment can undergo Secondary Fragmentation, producing a close low-mass binary, provided 
only that the primary fragment is spinning. Secondary Fragmentation is therefore most likely in primary fragments formed in 
the outer parts of circumstellar discs (since such fragments inevitably spin), and this could explain why a brown dwarf in a wide 
orbit about a Sun-like star has a greater likelihood of having a brown-dwarf companion than a brown dwarf in the field - as 
seems to be observed. Moreover, we show that binary brown dwarfs formed in this way can sometimes be ejected into the field 
without breaking up. 
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condensed out into seriously dense objects with masses in the 
ranges < 0.01 and/or > 100 M^, but - as far as we can judge 
- they have not. It is therefore appropriate to ask why gravi- 
tational condensation, under a wide variety of circumstances, 
selects objects in the mass range 0.01 to 100 M^. 

1.1. The maximum mass for star formation 



The upper mass limit is usually attributed to radiation pressure. 
As a forming star accumulates matter, its luminosity grows 
faster than its mass, and eventually - around 10 M^, in a spher- 
ically symmetric model - the outward force of radiation pres- 
sure becomes greater than the inward force of gravity, making 
further accretion difficult, but not impossible. More massive 



1. Introduction 

The structure, physical properties and appearance of the 
Universe are strongly influenced by the fact that most of the 
seriously dense baryonic condensations in it are objects with 
masses in the range 0.01 to 100 M^^, i.e. stars. By seriously 
dense we mean 



2 [47reJ 



^ 1 g cm 



(1) 



where is the critical density separating thermally ionised 
plasma from pressure-ionised plasma. The Universe would be 
very different if comparable amounts of baryonic mass had 
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stars can form, but their formation requires increasingly ex- 
treme and contrived circumstances. One possibility is that the 
inflow has already acquired such a high ram-pressure before it 
encounters radiation pressure that it cannot be reversed (McKee 
& Tan 2003). Alternatively, if collapse starts from sufficiently 
dense and homogeneous conditions, the accreting material is 
already within the dust-destruction radius when the central star 
switches on, and so the accreting material does not experience 
the outward force of radiation pressure acting on dust (Edgar 
& Clarke 2003). Another possibility is that accretion is highly 
aspherical (for example, by virtue of being channeled through 
a disc), and radiation pressure is then released along the direc- 
tions not occupied by the inflow (for example, along the rota- 
tion axis of the disc); this is sometimes called The Flashlight 
Effect (Yorke & Sonnhalter 2002). The final possibility is that 
stars above ~ lOM^ form by merging of lower-mass stars in 
the extremely dense and short-lived central cores of forming 
star-clusters (Bonnell, Bate & Zinnecker 1998). 



1.2. The minimum mass for star formation 

The lower mass limit is usually attributed to the thermody- 
namics of interstellar matter (Hoyle 1953; Rees 1976; Low & 
Lynden-Bell 1976; Silk 1977; Boss 1988). In order for a con- 
densation to fragment, the Jeans mass in the condensation must 
continue to decrease, and therefore the gas in the condensation 
must remain approximately isothermal. This in turn requires 
that the condensation must radiate away the PdV work done by 
compression, as fast as it is generated. During the early stages 
of contraction the rate at which PdV work is done on the con- 
densation (i.e. the compressional heating rate 'H) is of order the 
thermal energy divided by the freefall time. The radiative cool- 
ing rate (X.) depends on the emissivity and the optical depth 
of the condensation; it may also depend on how well the gas 
is thermally coupled to the dust, and on the ambient radiation 
field. The first concern of this paper is with evaluating the con- 
ditions under which X ^ 'H, and hence the minimum mass for 
star formation by Primary Fragmentation. 



1.3. The formation of binary systems 

A question closely related to the origin of stars is the origin 
of their binary properties. The subset of brown dwarfs that are 
in multiple systems appears to be concentrated in two distinct 
regions of the parameter-space of multiple systems. If the pri- 
mary is a Sun-like star, the semi-major axis is usually large, 
^ 100 AU. If the primary is a brown dwarf, the semi-major 
axis is usually small, ^ 20 AU. Moreover, as pointed out by 
Burgasser, Kirkpatrick & Lowrance (2005), there seems to be 
a significant population of hybrid triple or higher-order multi- 
ple systems in which a close ($ 20 AU) brown dwarf pair is 
in a wide orbit (^ 100 AU) around a Sun-like star. The sec- 
ond concern of this paper is to explore a possible explanation 
for these binary statistics in terms of Secondary Fragmentation, 
i.e. fragmentation facilitated by dissociation. 



1.4. Plan of the paper 

In Section 12 we introduce the constitutive physics which we 
will use. In Section|3] we analyse hierarchical fragmentation of 
a three-dimensional medium, and estimate the resulting min- 
imum mass. In Section |3 we re-evaluate the notion of hier- 
archical fragmentation, and argue that it probably not a good 
paradigm for star formation. In Section|5] we discuss one-shot 
fragmentation of a shock-compressed layer, formed by two col- 
liding streams. We derive the resulting minimum mass, and ar- 
gue that this is a much more realistic scenario for contemporary 
star formation in turbulent molecular clouds. We also explain 
why one cannot assume a priori that a minimum-mass fragment 
is marginally optically thick. In Section|6]we discuss fragmen- 
tation of a circumstellar disc and argue that this too may be an 
important source of low-mass stars, but that it can only be ef- 
fective at large distances {R ^ 100 AU) from the central star. 
This has important implications for the genesis of exoplanets 
and brown dwarfs, and for the truncation of massive protostel- 
lar discs. In Section0we discuss the nonlinear thermodynamic 
effects associated with H, dissociation and the resulting Second 
Collapse of a protostar. We explain how these effects may result 
in the formation of close, low-mass binaries, with separations 
in good agreement with observation. We also show that close, 
low-mass binaries could be formed in this way in the outer parts 
of circumstellar discs and then ejected into the field. In Section 
|S]we summarise our main conclusions. 

1.5. A note on nomenclature 

We shall assume that, when star-forming gas is approximately 
isothermal (say d{n[T]/d{n[p] $ 0.1), there are always density 
fluctuations on a wide range of scales. A density enhancement 
which might develop into a star is in the first instance termed 
a proto-fragment. If a proto-fragment is gravitationally unsta- 
ble, and its contraction time-scale is shorter than the timescales 
of competing processes, it becomes a fragment. If a fragment 
can radiate sufficiently fast to stay approximately isothermal it 
becomes a prestellar condensation, and this condition defines 
the minimum mass for Primary Fragmentation. The implication 
is that a prestellar condensation is significantly denser (say, at 
least one hundred times denser) than the initial proto-fragment. 
A prestellar condensation can subsequently undergo Secondary 
Fragmentation during The Second Collapse (i.e. when dis- 
sociates). 

2. Constitutive physics 

2. 1 . Equation of state 

We will be mainly concerned with contemporary star forma- 
tion in the disc of the Milky Way, and therefore with gas whose 
chemical composition (by mass) is 70% molecular hydrogen, 
28% atomic helium, and 2% heavy elements (in the form of 
molecules - like CO - and dust). The mean gas-particle mass 
is therefore m ^ 4.0 x lO^^^'g, the isothermal sound speed is 
a = [kjlfh]^!^ ^ 0.06 km s"' [r/K]'^^ and flie number den- 
sity of H2 is ^ p/[4.8xl0"^'*g]. At gas-kinetic temperatures 
T % 100 K, the rotational levels of are not strongly excited, 
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and therefore the adiabatic exponent is y ^ 5/3. At higher tem- 
peratures (300 K ^ r ^ 30, 000 K) y and fh decrease due 
to the excitation of the rotational and vibrational levels of H, , 
the dissociation of H,, and finally the ionisation of H°; we will 
consider one of the consequences of these changes in Section 

m 

2.2. Blackbody fluxes 

For algebraic convenience, we will frequently interchange T, m 
and a - [k^T /niY^^. In particular, it will sometimes be useful 
to express the flux from a blackbody surface as 



F = o- 



15c2/23 



(2) 



2.3. Rosseland- and Planck-mean opacities 

We will also assume that the Rosseland- and Planck-mean 
opacities due to dust are - to order of magnitude - the same, 
and given by 



K 



(3) 



with a:, = 10 ■'cm^g ' and emissivity index P - 2 for T % 
100 K (i.e. in the far-infrared and submillimetre wavebands). 

3. Hierarchical Primary Fragmentation of a 
three-dimensional medium 

3.1. The Jeans Mass in a 3-D medium 

In a uniform 3-D medium, an approximately spherical proto- 
fragment of mass m, will only condense out if it is sufficiently 
massive, 

1/2 

(4) 



m, > m„ 



36 G^p 



or equivalently if it is sufficiently small and dense. 



r < r 



3 G m. 



36 G^ m2 



(5) 
(6) 



Here the subscript '3' records that we are considering a proto- 
fragment trying to condense out of a 3-D medium. 

3.2. The condensation time-scale for Primary 
Fragmentation of a 3-D medium 

The timescale on which a Jeans-unstable proto-fragment con- 
denses out of a 3-D medium is given approximately by 

, nl/2 n -1/2 



^COND 3 (^3 ) 



32Gp 



1 



2/3 



i.e. proto-fragments with mass » 'Wjea^ss condense out on 
a freefall timescale, whereas less massive proto-fragments take 
longer, and proto-fragments with mass m, < 'Wje^ns? ^'^^ 
ever. 



3.3. Hierarchical Primary Fragmentation of a 3-D 
medium 

The molecular-cloud gas from which stars are forming today 
in the Milky Way is expected to be approximately isothermal, 
with T ~ lOK, as long as it can radiate efficiently via molec- 
ular lines and/or dust continuum emission. Therefore it has 
been argued, following Hoyle (1953), that star formation pro- 
ceeds in molecular clouds by a process of hierarchical Primary 
Fragmentation. An initially low-density massive cloud - which 
is destined to form a proto-cluster of stars - satisfies condition 
(0} and starts to contract. Once its density has increased by a 
factor the Jeans mass m^^^^^^^ is reduced by a factor 
and hence parts of the cloud can condense out independently, 
thereby breaking up the cloud into ^ / sub-clouds. Moreover, 
as long as the gas remains approximately isothermal (strictly 
speaking, as long as a remains approximately constant), the 
process can repeat itself recursively, breaking the cloud up into 
ever smaller 'sub-sub. ..sub-clouds'. 

3.4. The Opacity Limit for hierarchical Primary 
Fragmentation of a 3-D medium 

The process ends when the smallest sub-sub... sub-clouds are 
so optically thick, and/or are collapsing so fast, that the P dV 
work being done on them cannot be radiated away fast enough, 
and they start to heat up. It is normally presumed that this de- 
termines the minimum mass for star formation (e.g. Rees 1976; 
Low and Lynden-BeU 1976; Silk 1977; Boss 1988), and the ter- 
mination of hierarchical Primary Fragmentation in this way has 
traditionally been referred to as The Opacity Limit. Masunaga 
and Inutsuka (1999) point out that the gas does not have to be 
optically thick at The Opacity Limit. 

3.5. The compressional heating rate for a spherical 
fragment 

To estimate the minimum mass for hierarchical Primary 
Fragmentation of a 3-D medium, m^,,^,, we first formulate the 
P dV heating rate for a spherical fragment, neglecting the back- 
ground radiation field. 



K = -P. 



dV, 
dt 



3m^a^ dr^ 
r, dt 



3m.a 



Gm, 



1/2 



(8) 



Here we have obtained the final expression by assuming that 
the collapse is dynamical and putting drjdt [Gmjr^y^^. 

3.6. The radiative cooling rate for a spherical fragment 

The maximum radiative luminosity of a spherical fragment is 



(7) £, - 



3 [f^(T) + f-HT)Y 



where 



3 K„(T) 



(9) 



(10) 
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are the Rosseland- and Planck-mean optical depths. We empha- 
sise that T is the mean internal temperature of the fragment, 
and not necessarily the surface temperature. The treatment of 
radiation transport in Eqn. (|5Jl is based on the asymptotic forms 
for radiative diffusion in the optically thick limit (f^{T) » 1), 
and for local emissivity in the optically thin limit (fp(r) <s: 1). 
It is justified in more detail in Appendix A. Eqn. (|9jl gives a 
maximum net luminosity because it neglects the background 
radiation field (see below). 

3. 7. The Opacity Limit for hierarchicai Primary 

Fragmentation of a 3-D medium, assuming f ~ 1 

If we follow Rees (1976) and assume that the fragment is 
marginally optically thick, we can put [f^(r) + '(r)j ^ 2, 
and then - using Eqn. (|3 - the requirement that X, ^ "H, re- 
duces to 

3652 Gm^c^/i*']'^^ 



28;ri2 



311 5« Gh 



1/7 



(11) 



(12) 



This will give a conservative minimum mass, because, in as- 
suming that the fragment is marginally optically thick, we are 
maximising its luminosity - all other things (m^, and T) be- 
ing equal. 

3.8. Tine minimum mass for hierarclnicai Primary 
Fragmentation of a 3-D medium, assuming f ~ 1 

Conditions (|5} and (lll> require r^^ 



■JEANS, > '•cool, and hence 



283 



PLANCK 



1/2 



(13) 



Here «Vanck - [hc/G]^^^ = 5.5 x 10"^ g is the Planck mass, 
and m cannot differ greatly from (the proton mass), so m^^^^ 
is essentially the Chandrasekhar mass (~ irr' Inp-) times a 

^ PLANCK ' p ^ 

factor [a/c]'^^ ~ 10"^. We note the relatively weak dependence 
of m^^^^ on T (oc r^/^) and the relatively strong dependence on 

in (oc m"^^^). 

For contemporary local star formation, we substitute in ^ 
4.0xl0"2'*gandfl ^ 1.8 xlO"* cms (corresponding to molec- 



ulai- gas at r ^ 10 K) and obtain 



0.004 M„. 



3.9. Allowing for radiation transport effects in 
hierarchical Primary Fragmentation of a 3-D 
medium 

There is no a priori reason why the limiting fragment should be 
marginally optically thick (Masunaga & Inutsuka 1999), and 
therefore we cannot necessarily put [f^(r) + Tp"'(r)] - 2. 

If we assume that the fragment is very optically thin, we 
must put [T^(r) + Tp"'(r)] - fp '(r), and then after some te- 
dious but straightforward algebra we obtain 



31/25 



23 7r2 Gai'^a^-KjT) 



(14) 



Conversely, if we assume that the fragment is very optically 
thick, we must put [fp(r) + f~^{T)^ ^ ^^T), and then after 
some more tedious but straightforward algebra we obtain 



5;r* c^h^a^K„iT) 



1/3 



2733/2 G5w4 ■ (^^^ 

In general one should use both Eqn. (I14> and Eqn. jl5> . and 
then adopt whichever gives the larger value of /«M,t^3. 

In the case of contemporary local star formation with T ^ 
10 K, in ^ 4.0x10-^'^ g, a ^ l.SxlO^kms-', /f, ^ lO^^cm^g-' 
and /3 - 2, both expressions (Eqns. and I15> again give 
m„,^, ~ 0.004 Mj,, but this is because - coincidentally - frag- 
ments of this mass are marginally optically thick (f ^ 1) at 
their inception. 

4. Problems with hierarchical Primary 
Fragmentation of a 3-D medium 

The above analysis of three-dimensional hierarchical Primary 
Fragmentation overlooks several effects which all act to in- 
crease m„,fj3. Indeed, three-dimensional hierarchical Primary 
Fragmentation probably does not work. There is no conclusive 
evidence that hierarchical Primary Fragmentation happens in 
nature, nor does it seem to occur in numerical simulations of 
star formation. There are three main reasons. 

4.1. Merging 

First, Eqn. clearly shows that at each level of the hierarchy, 
a proto-fragment condenses out more slowly than the larger 
parent fragment of which it is part, by virtue of the fact that 
the proto-fragment is always less Jeans unstable than the par- 
ent fragment. Therefore a proto-fragment is unlikely to survive 
as a distinct entity. Instead the overall contraction of the par- 
ent fragment will usually cause a proto-fragment to merge with 
neighbouring proto-fragments. 

4.2. Accretion 

Second, if there are no neighbouring proto-fragments with 
which to merge, then, even if an isolated proto-fragment starts 
off with a mass m, ~ nij^^^^,, it will subsequently grow by a 
large factor due to accretion. Therefore, before its contraction 
approaches freefall and it separates out from the background, 
its final mass will have become much greater than wJjeans, ■ 

Although the following estimate is only indicative, we 
adopt the formula for Bondi accretion (Bondi, 1952) . This 
gives a growth rate 



dm, 
dt 



e^''-KG^pm^ 



m. 



23 61/2 



■■JEANS 3 



22 



(16) 



where t^^ - [37r/32Gp] ' is the freefall time in the background 
medium of the parent fragment. The implication of Eqns. (0 
and (I16> is that a proto-fragment with w, ~ 'Wjeans 3 takes several 
freefall times to condense out, and during one freefall time it 
increases its mass by more than an order of magnitude. 
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4.3. Back-warming 

Third, individual proto-fragments will be back-warmed by the 
ambient radiation field from other cooling proto-fragments, 
which in principle fill a significant fraction of the celestial 
sphere, and again this will reduce their net luminosity X,, 
thereby inhibiting fragmentation. 



5. Two-dimensional one-shot Primary 

Fragmentation of a shoclt-compressed layer 

5.1. Colliding flows, and star formation in a crossing 
time 

In reality, 3-D hierarchical Primary Fragmentation may be an 
inappropriate paradigm for star formation in molecular clouds. 
There is growing evidence that star formation in molecular 
clouds proceeds very rapidly, essentially 'in a crossing time' 
(Elmegreen 2000), i.e. whilst the molecular cloud is being as- 
sembled. In this scenario, star formation occurs in molecular 
clouds where two turbulent flows of sufficient density collide 
with sufficient speed - and therefore sufficient ram pressure - 
to produce a gravitationally unstable shock-compressed layer 
which then fragments to produce prestellar condensations. This 
mode of Primary Fragmentation is two-dimensional because 
the motions which assemble a prestellar condensations out of 
a shock-compressed layer are, initially, largely in the plane of 
the layer. It is one-shot in the sense of not being hierarchical. 

5.2. Geometric considerations 

Strictly speaking, our discussion of shock-compression should 
not be limited to layers, but should also consider other geome- 
tries, in particular filaments and isolated globules. However, we 
believe that the Primary Fragmentation of a shock-compressed 
layer captures the essential elements characterising turbulent 
fragmentation, viz. the dynamic interplay between ram pres- 
sure and self gravity in producing prestellar condensations 
from interstellar gas which is able to keep cool radiatively. 

It might be argued that there is limited observational evi- 
dence for shock compressed layers, but this is probably a se- 
lection effect. Such layers will only be easy to discern under 
the relatively rare circumstance that they are viewed edge-on 
— and then they will be hard to distinguish from filaments. 

Moreover, when a shock-compressed layer becomes grav- 
itationally unstable, it fragments first into filaments, and then 
into individual prestellar condensations (e.g. Turner et al. 1995, 
Whitworth et al. 1995, Bhattal et al. 1998). The separation be- 
tween adjacent filaments in the layer is then essentially the 
same as the separation between prestellar condensations in 
the same filament (and is given by the fragmentation length, 
2'"fragt defined in Ean.Bllbelow). This is a generic property of 
the gravitational fragmentation of a flattened structure, and is 
also seen in cosmological simulations. 

Filaments can also be formed directly if more than two tur- 
bulent flows collide more-or-less simultaneously, and this is 
seen in simulations of interstellar turbulence (e.g. Hennebelle 
& Audit 2005; Vazquez-Semadeni et al. 2006). A filament 



formed in this way may then be dense enough to fragment grav- 
itationally into a chain of prestellar condensations. 

More rarely, convergent turbulent flows may be sufficiently 
well focussed to form a single compressed prestellar condensa- 
tion which then collapses in isolation (Whitworth et al. 1996). 

However, we will consider here only the case of a shock 
compressed layer, because it is the simplest case to treat ana- 
lytically, and probably also the most common case. 

5.3. Growth of a shock-compressed layer 

Specifically, we limit our discussion to the generic case of two 
flows, both with uniform pre-shock density p, colliding head- 
on at relative speed v. We assume that the effective isothermal 
sound speed, a, in the post-shock gas of the layer satisfies a «; 
v/2, or in other words, that the Mach Number of the accretion 
shock bounding the layer is large: 



M - - » 1 

Za 



(17) 



Therefore the surface-density £ of the layer grows according to 

nt) ^ pvt- (18) 

its density is ~ pAf' ~ and its half-thickness is 

2fl2f 

z(t) - . (19) 



5.4. Structure of a shock-compressed layer 

The layer is contained by the ram pressure of the inflowing gas, 
rather than its self-gravity, as long as ^ pv^/4 » G2^(f), 
which reduces to 



1 



2[Gp] 



1/2 ■ 



(20) 



Moreover, there is plenty of time for the layer to relax to 
hydrostatic equilibrium, as long as its sound crossing time, 
fsc(0 - z(t)la, obeys the condition t^^it) <k t. This condition 
reduces to v » 2a, and is therefore automatically satisfied (see 
EQn.ll7>. Thus, until the layer fragments, it has a rather flat den- 
sity profile, i.e. the density throughout the layer is ~ pv^/4a^. 

5.5. How to estimate the scale on which a 

shock-compressed layer undergoes Primary 
Fragmentation 

Before proceeding to estimate the minimum mass for this case 
('"mini)' is appropriate to sketch the steps which our analy- 
sis will take. Because the surface-density of the layer increases 
monotonically with time (Eqn. [TSt. the Jeans mass m^^^^^^it) 
(i.e. the minimum mass for a proto-fragment which could con- 
dense out of the layer, given sufficient time) decreases mono- 
tonically with time (Eqn. I23> . Furthermore, it turns out that 
at time t proto-fragments with mass m^^^^^^,(t) ^ 4OTjg_^^s2(0 
condense out fastest of all (i.e. faster than both smaller and 
larger proto-fragments), and on a time-scale fFASTEST2(0 (Eqn. 
l28l . Therefore we hypothesise that non-linear fragmentation 
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of the layer takes place at time t-^f^^-. when the condensation 
time-scale of the fastest condensing proto-fragment is of or- 
der the elapsed time, i.e. t^^,^,,,{t^^^,) ^ t^^, (Eqn. H^, 
and this determines the mass of a typical fragment, m^^^^^ ^ 
'"fastest2(Wg2) (Eqn. 1^ . We note that this hypothesis has 
been corroborated by numerical simulations using a variety of 
numerical codes (Whitworth et al. 1995; Turner et al. 1995; 
Bhattal et al. 1998; Dale, Bonnell & Whitworth 2006). Finally, 
is the smallest value of m^^^, for which the fragment can 
radiate sufficiently fast to keep cool as it condenses out. 



5.6. The Jeans mass in a shock-compressed layer 

A disc-like proto-fragment of radius r^, in a shock-compressed 
layer, evolves according to 



dV, 
df^ 



-7rGS(0 + —. 



(21) 



where the first term on the righthand side of Eqn. (12 1> repre- 
sents self-gravity, and the second term represents hydrostatic 
support. It follows that at time t only proto-fragments with 



r ^ r (t) ^ 



nGY.{t) nGpvt ' 



(22) 



2 JEANS 2 V ^ JEANS 2 ^ ^ ^ ^ 



7rG2S(0 nG'-pvt ' 



(23) 



can start to condense out. 

Here the subscript '2' records the fact that we are consider- 
ing a proto-fragment trying to condense out of a 2-D layer. 

5.7. The condensation time-scale for a proto-fragment 
in a shock-compressed layer 

It is important to note that, in a shock-compressed layer, the 
Jeans mass is a function of time; it decreases monotonically 
as the surface density of the layer increases. We now need 
to evaluate the condensation time-scale for a Jeans-unstable 
proto-fragment, in order to determine the stage at which such a 
proto-fragment actually has sufficient time to condense out of 
the growing layer. 

At time t, the time-scale on which an unstable proto- 
fragment condenses out is given by 



2('-2>0 



1/2 



7rGS(f) a' 



2^-U2 



nGpvt 



At) 



At) 



-1/2 



(24) 



or, in terms of m. 



fcOND2('«2'0 



nGpvt 



1/2 



'«JEANS2(0 



-1/2 



(25) 



5.8. The fastest condensing proto-fragment in a 
shock-compressed layer 

It follows that, at time t, the fastest condensing proto-fragment 
has 



2(0 ^ 2r„_,(f) 



2fl2 



nGpvt 



(26) 



2 

'"fASTEST2(0 ^ ^^^5X2^*^^) ^ ^OT^^ANSjCO ^ j^Ql^^f ' "^^T) 

and that the fastest condensing proto-fragment condenses out 
on a time-scale 



2(0 



2a 



nGpvt 



(28) 



We note that it is a generic property of the Primary 
Fragmentation of a layer that there is a finite size of proto- 
fragment which condenses out faster than both smaller proto- 
fragments and larger proto-fragments (e.g. Larson 1985). 

5.9. The epoch of non-linear Primary Fragmentation 
for a shock-compressed layer 

^From Eqn. J28> it follows that the accumulating layer will start 
non-linear fragmentation once fpASTESTTCO ^ t, i.e. at 



t ^ t„ 



2 a 



nGpv 



1/2 



(29) 



5.10. The scale of non-linear Primary Fragmentation 
for a shock-compressed layer 

Once non-linear Primary Fragmentation occurs in a shock- 
compressed layer, the characteristic scale of the fragments 
which form is 



= z, 



■FASTEST 2 '- FRAG 



(^FRAG^) 
2 ) 

2 ) ~ 



8fl^ 



nG^ pv 



1/2 



2a^ 



nGpv 
8a^ 



1/2 



1/2 



nGpv^ 



(30) 
(31) 
(32) 



5.11. Distinctive features of the Primary 

Fragmentation of a shock-compressed layer 

We assume that non-linear Primary Fragmentation occurs as 
soon as it is possible, i.e. at f ~ ^frag-' ^^'^ results in frag- 
ments of mass ~ ni^^^^. This presupposes that there is enough 
substructure in the pre-shock gas on scales ~ '"fragt 
seed the fastest condensing proto-fragments. With this assump- 
tion, we note the following distinctive aspects of the Primary 
Fragmentation of a shock-compressed layer. 

(i) The characteristic mass, mp^^^,,, and initial radius, rp^^G" 
of a fragment both decrease monotonically with the mass flux, 
pv, in the colliding flows. 
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(ii) Fragments are flattened at their inception, with an as- 
pect ratio approximately equal to the Mach number of the ac- 
cretion shock bounding the layer: 



:^ ^ M » 

2a 



1 



(33) 



(iii) Wp^^j,, is not simply the 3-D Jeans mass {m^^^^^^; Eqn. 
13 evaluated at the post-shock density (~ pv^/4a^). It is larger 
by a factor - 12M^^^/7t\ 

(iv) The shock-compressed layer is still contained by the 
ram pressure of the inflowing gas, rather than by self-gravity, 
when its non-linear Primary Fragmentation starts, provided that 
tmAoi < l/2[Gp]'^^ (see Eqn.|20j, i.e. provided that M > 4/;r, 
which is satisfied automatically (see Ean.ll7>. 

(v) The assumption that matter is still flowing into the ac- 
cretion shock bounding the layer at time fpn^Q, requires that the 
linear extent, L, of the pre-shock gas normal to the shock ex- 
ceed t^^^,v/2, i.e. 



T > T ~ 



2tiG p 



1/2 



(34) 



(vi) If the pre- and post-shock gases both subscribe to a 
Larson-type scaling relation, so that the effective sound speed, 
QgPP (including the contribution from the turbulent velocity dis- 
persion), scales with the mean density, p, according to oc 
p then the pre-shock gas can have the required extent, L, 
and still be gravitationally stable, provided 



[Gp] 



1/2 



4Gp 



1/2 



(35) 



(vii) Since L^^^/L„^ ^ [nM/2y''^ > 1, there is a range 
of L-values for which the pre-shock gas is gravitationally sta- 
ble, but the shock-compressed layer still undergoes Primary 
Fragmentation whilst it is accumulating. 

(viii) As in hierarchical Primary Fragmentation, a fragment 
in a shock-compressed layer will only condense out if it is able 
to remain approximately isothermal by radiating efficiently. 



5. 12. The compressional heating rate for a fragment in 
a shock-compressed layer 

The compressional heating rate for a flattened disc-like frag- 
ment in a layer is 



■H, = -P. 



pv^ 2 7Tr^Z2 2a^ 



dt 4 f^o^o, G 



(36) 



where we have obtained the final expression by substituting 

^2 ~ ^FRAG2' ^2 ~ ^FRAG2' '^^^ ^COND2 ~ ^FRAG2' 

5. 13. The radiative cooling rate for a fragment in a 
shock-compressed layer 

The radiative cooling rate for a flattened disc-like fragment is 



[f,(T) + f-HT)] 



(37) 



where the optical depths are now given by 



(38) 



A detailed justification for Eqns. ( I37> and ( I38> is given in 
Appendix A. 

5.14. The Opacity Limit and the minimum mass for 
Primary Fragmentation of a shock-compressed 
layer 

The requirement that X, ^ "H, then reduces to a limit on the 
mass flux, pv, in the colliding flows. 



2 (,3 



pv 



[T,(r) + f-i(r)] 



(39) 



If the fragment is marginally optically thick, we can set 
[f,(r) + Vi(r)]-2,andthen 



171 ^ lU 



[30] 



1/2 



1/2 



If the fragment is optically thin, we obtain 



15 



2 1.3 



C^h 



47r5 Gm'^a^K^iT) ' 
If the fragment is optically thick, we obtain 



60 c^h^a^KjT) 



71 



G5 



nr 



1/3 



(40) 



(41) 



(42) 



In general one should evaluate both Eqn. J41l i and Eqn. 
(142 > ■ and then use whichever gives the larger value of m^^^^. 

In the case of contemporary local star formation with T ^ 
10K,m ^4.0x10-24 ^ 1.8x10'* kms-i,/^, ^ lO^^cm^g-i 
and/? ^ 2, (Eqns.|40|il]and|43 all give m„,^, ~ 0.001 M^, but 
this is because once again the minimum fragment is - coinci- 
dentally - marginally optically thick. 

We stress that this really is a coincidence. The minimum- 
mass fragment does not have to be marginally optically thick. 
At any temperature, there is a critical opacity. 



15 
(2^ 



G2c4/i« 



1/4 



(43) 



Under circumstances where the actual opacity falls below the 
critical value, k^{T) < k^^^^(T), the minimum-mass fragment is 
optically thin (and cannot cool fast enough because it is opti- 
cally thin). Conversely, when K^iT) > ^^^.^(r), the minimum- 
mass fragment is optically thick (and cannot cool fast enough 
because it is optically thick); parenthetically, in this regime, 
somewhat larger fragments which can cool fast enough are also 
optically thick, but still (approximately) isothermal. 

For r ^ lOK and OT ^ 4 X 10"^"* g, we have k^^„(T = 
10 K) ^ 0.11 cm^ g"', which just happens to be close to our es- 
timate of the opacity in local star forming gas, K^^iT - 10 K) ^ 
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0.10 cm- g"' . On the assumption that the amount of dust is pro- 
portional to the metallicity, Z, but the intrinsic mix of grain 
properties is universal, we have 



^crit(-^) 



0.02 



10 K 



15/4 



4 X 10-24 g 



1/4 



(44) 



Thus, in regions like the outer parts of the Galaxy, where the 
metallicity and/or the temperature are lower, the optically thin 
expression should be valid. Conversely, in regions like the inner 
parts of the Galaxy, where the metallicity and/or the tempera- 
ture are higher, the optically thick expression should be valid. 

We note also that the minimum mass for two-dimensional, 
one-shot Primary Fragmentation, m^^^^ (Eqns. ^] |^ & I42> . 
depends on physical variables in exactly the same way as 
the minimum mass for three-dimensional, hierarchical Primary 
Fragmentation, (Eans. [T3lll4l &[T5l. This too appears to 
be somewhat fortuitous, given that the expression for m^g^^j, 
(Eqn. 0} is quite diferent from that for '«jeans2 (Eqn. I23> or 
'"frag: (Eqn.|30j. 

5.15. Advantages of two-dimensional, one-shot 
Primary Fragmentation of a layer 

In addition to delivering a smaller minimum mass, two- 
dimensional one-shot Primary Fragmentation bypasses the 
three problems associated with three-dimensional hierarchical 
Primary Fragmentation which we discussed in Section^] 

First, there is little likelihood of neighbouring fragments 
merging, since fragments with mass ~ '«frag2 condense out of 
the layer faster than larger structures (see Eqn. i25\ and Larson 
(1985)). 

Second, accretion is less problematic. Boyd & Whitworth 
(2004) have analysed in greater detail the radiative cooling of a 
fragmenting layer, using a two-dimensional model, and taking 
into account both on-going accretion (as matter continues to 
flow into the layer) and the energy dissipated in the accretion 
shock. They find that the rate at which energy is dissipated in 
the accretion shock is comparable with the rate at which the 
gas in the fragment is heated by compression. For T ^ 10 K, 
the minimum mass is ~ 0.0027 M^. This fragment starts with 
mass ~ 0.0011 M^, but continues to grow by accretion as it 
condenses out. The minimum mass can be reduced further by 
decreasing T. 

Third, there is little backwarming, because there are no 
other fragments filling the part of the celestial sphere into 
which a fragment radiates (i.e. perpendicular to the shock- 
compressed layer). 

6. Primary Fragmentation of a circumstellar disc 

Another scenario which may be of more relevance to con- 
temporary star formation than three-dimensional hierarchical 
Primary Fragmentation is the Primary Fragmentation of a mas- 
sive circumstellar disc. Massive disc-like structures are ob- 
served around some Class and Class I protostars (e.g. Eisner 
et al. 2005; Eisner & Carpenter 2006). However, they are quite 
rare, and we therefore infer, either that discs seldom form, or 



that they are short-lived. Numerical simulations support the lat- 
ter inference. For example, in simulations of the collapse and 
fragmentation of intermediate- and low-mass turbulent cores 
(Bate et al. 2003; Goodwin et al. 2004), massive disc-like struc- 
tures form quite commonly around the initial (primary) proto- 
stars, but then fragment to form secondary protostars before 
they can relax to equilibrium. Although massive circumstellar 
discs may be transient structures, we can only analyse their sta- 
bility analytically if we assume that they are sufficiently regular 
to be approximately azimuthally symmetric. 

6.1. The Toomre criterion for Primary Fragmentation 
of an equilibrium circumstellar disc 

A necessary, but not sufficient, condition for the Primary 
Fragmentation of an azimuthally symmetric, equilibrium cir- 
cumstellar disc is that the surface density S be sufficiently large, 

a(R) e(R) 



I.{R) > l^^^^^^iR) ^ 



7tG 



(45) 



(Toomre 1964). Here R is distance from the primary protostar at 
the centre of the circumstellar disc, a{R) is the local isothermal 
sound speed. 



e(R) 



1r 



R- 



+ 4Q.^(R) 



1/2 



(46) 



is the local epicyclic frequency', and Q.(R) is the orbital angular 
speed. 

6.2. Criterion for Primary Fragmentation of a 
non-equilibrium circumstellar disc 

However, since in simulations of star formation circumstel- 
lar discs fragment before they have time to approach dynam- 
ical equilibrium, it is useful to have a criterion for Primary 
Fragmentation which can also be applied to non-equilibrium 
circumstellar discs, where the epicyclic frequency is not prop- 
erly defined. Such a criterion can be derived, based on the lo- 
cal vorticity, a){R), which is defined even for non-equilibrium 
discs; the only constraint is that the vorticity must be slowly 
varying on the scale of a proto-fragment, so that it gives a re- 
alistic measure of the spin of the proto-fragment. We note that, 
in the case of a Keplerian disc, the new constraint (Eqn. 
below) is exactly equivalent to the old one, because in this cir- 
cumstance 2co(R) - Q.{R) - eiR). Therefore the results we 
present subsequently are not affected by which criterion we use 
(Eqn.|45lor|lH}. 

Consider a small disc-like proto-fragment of radius with 
its centre at distance R from the centre of the circumstellar disc. 
Assume that the extent of the disc-like proto-fragment, ~ 2rQ, 
is much smaller than the extent of the circumstellar disc of 
which it is part, ^ 2R. Radial excursions of the proto-fragment 
evolve according to 



dV, 
df2 



- 7tG'E(R) 



a\R) 



u\Ry^ . 



(47) 



' We have broken with convention in calling the epicyclic frequency 
e, rather than k, simply to avoid confusion with the opacity. 
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In Eqn. J47l i. the first two terms on the righthand side repre- 
sent, respectively, self-gravity and hydrostatic support (just as 
in Ean.l21>. The third term on the righthand side of Eqn. j47t 
represents centrifugal support. The subscript 'D' records the 
fact that we are considering a proto-fragment trying to con- 
dense out of a circumstellar disc. 

The righthand side of Eqn. ( I47> has its minimum when r-^ ^ 
a{R)/cL>(R). This minimum is negative - signifying instability 
against contraction - only if 



2 aiR) aj(R) 

7g 



(48) 



We re-iterate that this criterion for Primary Fragmentation of 
a general disc (Eqn. I48t reduces to The Toomre Criterion 
(Eqn. 145 > if we substitute 2lo{R) = e(R), as appropriate for a 
Keplerian disc. 

6.3. The condensation time-scale for a proto-fragment 
in a circumstellar disc 

^From Eqn. ( I47> . the condensation time-scale for a proto- 
fragment in a circumstellar disc is given by 

1/2 



7tGI.(R) a\R) 2 



co\R) 



-1/2 



(49) 



6.4. The fastest condensing proto-fragment in a 
circumstellar disc 

^conddC-'^' '"d) has a minimum for ^ 2a^(R)/nGl.(R), and 
therefore the radius, mass and growth time-scale of the fastest 
condensing proto-fragment at radius R are given by 



'Wfastestd(^) ^ 4m_ 



2a\R) 
7rGI.{R) ' 
4-a\R) 



^TESTD 

(R) ~- 

If we define 

2 a{R) aj(R) 




Q(R) = 



nClLiR) 



(50) 
(51) 

(52) 
(53) 



the condition for Primary Fragmentation of a general disc (Eqn 
BSl becomes Q{R) < 1, and Eqns. (|50j to (|52} reduce to 



^FASTEST D 

{R) - 

^FASTEST D (^) ~ 
FASTEST D 

{R) - 



(jj(R) ' 
2Q(R)a\R) 



GcoiR) 



1 



IQ HR) - if^ oj(R) 



(54) 
(55) 
(56) 



6.5. The typical fragment mass in a dynamically 
forming circumstellar disc 

Eqn. ( 145 > is also approximately the condition for spiral modes 
to develop in an equilibrium circumstellar disc, and these will 
have the effect of redistributing angular momentum. As a re- 
sult, the inner parts of the circumstellar disc may simply ac- 
crete onto the central primary protostar, and the outer parts 
may simply disperse without fragmenting. The implication is 
that Primary Fragmentation is much more likely if - as in the 
numerical simulations of Bate et al. (2003), Goodwin et al. 
(2004), Hennebelle et al. (2004) - the approach to instability 
is dynamic, rather than quasistatic. 

To explore this dynamic situation, we assume that the cir- 
cumstellar disc grows in mass rather rapidly, and is thereby 
launched directly into non-linear fragmentation with Q ~ 0.5 
(rather than Q edging gradually downwards past unity, and in- 
stability having first to grow slowly through a linear phase). 
The radius, mass and condensation time-scale of a typical frag- 
ment are then 



^FRAGD (R) 
m^AaJR) 
^FRAG D (R) 



a(R) 
2MK) ' 

a^{R) 
Gu)(R) ' 
1 

yi^a){R) 



(57) 
(58) 
(59) 



6.6. The Gammie criterion for fragmentation of a 
circumstellar disc 

A second necessary, but not sufficient, condition for the 
Primary Fragmentation of a circumstellar disc is that fragments 
be able to cool radiatively on a dynamical time-scale (Gammie 
2001); this is just The Opacity Limit under a slightly differ- 
ent guise. Specifically, Gammie suggests that the cooling time- 
scale for a fragment, t^^^^, must satisfy t^^^^^ ^ t^^^^^l2, where 
fgi^g,^ is the local orbital period. This condition has been corrob- 
orated by numerical simulations (Rice et al. 2003). If a frag- 
ment cannot cool sufficiently fast, it is likely to bounce and be 
sheared apart (e.g. Cai et al., 2006). 

6.7. The compressional heating rate for a fragment in 
a circumstellar disc 

The compressional heating rate for a fragment at radius in a 
circumstellar disc is 



'H,{R) 



3m,^,,^(R)a\R) 3^'^ a\R) 



2 ^FRAG D (^) 



2G 



(60) 



6.8. The radiative cooling rate for a fragment in a 
circumstellar disc 

The radiative cooling rate for a fragment in a circumstellar disc 

is 

A-nr^ (R)o-,,T(R)'^ 

PRAr. n*- SB 



-C,iR) - 



[f,(T) + f-HT)] 



(61) 
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and the optical depths are given by 

'"fragd(^)'^m(7') 



temperature for Primary Fragmentation of a circumstellar disc: 



fJT) - fJT) 



(62) 



These are essentially the same expressions as were invoked for 
Primary Fragmentation of a shock-compressed layer (Eqns.lTTl 
andl38>, since here too we are dealing with a disc-like fragment 
which radiates mainly through its flat surfaces. See Appendix 
A for a detailed justification of Eqns. i6l\ and i62\ . 

6.9. The Opacity Limit for Primary Fragmentation of a 
circumstellar disc 

The requirement that X.^{R) <: 'H^iR) therefore reduces to 



(63) 



aHR) 35/2 5 c2 [f^(T) + f -i(r)] ' 

6.10. A specific circumstellar-disc model 

To make the discussion more specific we consider a quasi- 
Keplerian disc around a Sun-like star having mass M*, lumi- 
nosity L*, and hence 



7 



to(R) ~ 10"' s 



T(R) ~ 300 K 





1/2 


R 


-3/2 






IauJ 







1/4 


R 


-1/2 


K 




IauJ 





a(R) 



1 kms 





1/8 


R 


-1/4 






IauJ 





(64) 
(65) 
(66) 



The typical fragment at radius R then has radius, mass and 
optical depth 



r^,,,(R) ~ 0.03 AU 





-1-1/2 


L 


*■ 


1/8 


R 




5/4 


(67) 




























L 









AUJ 










-1/2 


L 


★ 


3/8 




R 


3/4 


Mo 




















, (68) 




Me 












LAUJ 




1/2 




5/8 


r R ^ 


-11/4 


(69) 
































LAUJ 







f,,,^,{R) ~ 2x10' 



6.11. The minimum mass for Primary Fragmentation 
of a circumstellar disc 

^From Eqn. ( I69> we can assume that the fragment is optically 
thick, with [fJT) + f;'(T)] - f^^^^^iR). Eqn. ^ then re- 
duces to the form 



R ^ R.. 



150 AU 



1/3 



(70) 



i.e. prestellar condensations can only form in the outer parts of 
a circumstellar disc, because only in the outer parts of a disc 
can such condensations radiate fast enough. 

Substituting from Eqn. ilOt in Eqns. ( I67t . ( I68> and i65\ . 
there is a minimum initial condensation radius, a minimum ini- 
tial condensation mass, and a maximum initial condensation 



Vnd ~ 16 AU 





-1/12 




1/8 


ImJ 









0.003 



^MAXD ~ 25K 





-1/4 




[mJ 






-1/6 




1/4 









3/8 



(71) 
(72) 
(73) 



Similar conclusions (Eqns. l7Ull71ll72l andl73t were reached 
by Rafikov (2005) using closely related arguments. 

6.12. The Brown Dwarf Desert 

Although we have considered a disc with specific surface- 
density, velocity and temperature profiles, the final result is not 
very sensitive to these assumptions, as evidenced by the rela- 
tively small exponents in Conditions J70> . J71> , J72> and (I73> . 
We can allow for some variance in the disc parameters by relax- 
ing these conditions slightly. In particular, we adjust Condition 
(|7U)to 



R ^ 100 AU. 



(74) 



If brown-dwarf companions to Sun-like stars are formed by 
Primary Fragmentation of discs, then the constraint that this 
can only occur at large radii, R ^ 100 AU, may explain The 
Brown-Dwarf Desert. 

6.13. Dissipation of massive extended protostellar 
discs 

Observational estimates of the specific angular momentum, 77, 
in star-forming cores (Bodenheimer, 1995; his Fig. 1) give 
77 10^' cm^ s"'. If deposited in orbit around a 1 star, this 
material should end up at radius 



R 



> 300 AU. 



(75) 



It is therefore noteworthy that massive circumstellar discs of 
this extent are quite rare. The simplest explanation is that 
they are short-lived, being converted rapidly into brown dwarfs 
and low-mass H-burning stars, on a dynamical timescale (~ 
lO'* years). 

6.14. Forming exoplanets by Primary Fragmentation 

Conversely, our analysis shows that Primary Fragmentation is 
strongly inhibited at small radii, say R ^ 30 AU, by the in- 
ability of a proto-fragment to radiate fast enough. Therefore it 
seems very unlikely that gas-giant exoplanets form frequently 
by Primary Fragmentation. This conclusion is in agreement 
with the analysis of Rafikov (2005), and with the simulations 
of Cai et al. (2006). In contrast. Boss (e.g. 2004) has suggested 
- again on the basis of numerical simulations - that gas-giant 
planets are able to form by Primary Fragmentation at small 
radii, because they can cool convectively at these radii. We cau- 
tion that Boss's interpretation of his numerical results is very 



A. P. Whitworth and D. Stamatellos: The minimum mass for star formation, and the origin of binary brown dwarfs 



11 



speculative. A proto-fragment condensing out gravitationally, 
on a dynamical timescale, does not have sufficient time to cool 
its interior by convection, because this would require the con- 
vective cells to migrate and disperse supersonically. What Boss 
may be seeing is a fragment undergoing an 'adiabatic bounce' 
prior to being sheared apart by differential rotation. 

7. dissociation and Secondary Fragmentation 

7.1. The spin angular momenta of prestellar 
condensations 

If low-mass prestellar condensations form by Primary 
Fragmentation in the outer parts of massive circumstellar discs, 
in the manner analysed in the preceding section, their subse- 
quent contraction is likely to be moderated by the rate at which 
they can lose spin angular momentum. This is because, in a 
marginally unstable disc, there is only a small range of unsta- 
ble fragment masses. These fragments are only just big enough 
for their self-gravity to overcome internal pressure, and only 
just small enough for their self-gravity to overcome centrifu- 
gal acceleration. Therefore such fragments are naturally born 
in a state where spin angular momentum makes a significant 
contribution to their subsequent dynamical evolution. 

Moreover, as a prestellar condensation contracts, and its 
moments of inertia decrease, its ability to loose spin angular 
momentum to the surroundings by gravitational torques de- 
creases, and the condensation becomes progressively more iso- 
lated - dynamically - from the rest of the disc. Recent MHD 
simulations of a protoplanet accreting from a protoplanetary 
disc (Machida, Inutsuka & Matsumoto, 2006) show that the 
protoplanet is able to lose angular momentum by launching 
bipolar jets. However, this loss of angular momentum does not 
occur on a dynamical timescale. Therefore we can still presume 
that condensations forming in discs by Primary Fragmentation 
derive significant support from rotation, and are quite flattened. 

A similar situation holds for low-mass prestellar conden- 
sations formed in other ways, for example by hierarchical 
Primary Fragmentation of a 3-D medium (Section |3j, or by 
one-shot 2-D Primary Fragmentation of a shock-compressed 
layer (Section|5}- All that is required for a condensation to be 
flattened by rotation is that its specific spin angular momentum, 
{, satisfy 



FRAG 



-1/6 



(76) 



We note that the prestellar condensation can have had a higher 
^- value when it first formed, but to reach density p without frag- 
menting, i must have been reduced to a value on the order given 
by Eqn. ( l76l . 

7.2. The Second Collapse 



When the temperature in a prestellar condensation reaches 
Jpigj, ~ 2, 000 K and the density reaches Pp,j,5 ~ 10"^ gem"-', 
first the vibrational degrees of freedom of start to be excited, 
and then starts to dissociate. As the temperature continues 
to increase, more than half the self-gravitational potential en- 
ergy released by contraction has to be invested in rotational 



and vibrational degrees of freedom, and then in dissociation of 
H,, rather than being invested in translational degrees of free- 
dom. As a result the pressure falls below the value required for 
quasistatic contraction and the condensation collapses (Larson 
1969). In the context of isolated protostars, this is normally re- 
ferred to as The Second Collapse. 

The Second Collapse is further promoted by the fact that 
these temperatures and densities correspond to The Opacity 
Gap, i.e. the regime where dust has sublimated and the H ion 
is not yet abundant, so there is only a very low residual opac- 
ity due to molecular lines. As a result, radiative energy leaks 
out of the condensation very rapidly, accelerating its collapse. 
However, this is a secondary factor when compared with H^ 
dissociation (see Whitworth et al. in preparation). 

We note that the critical value of the specific angular mo- 
mentum for H, dissociation is 



10 



19 2 - 

cm s 



0.1 M„ 



2/3 



(77) 



(where we have simply substituted p ~ p^^^^ in Eqn. I76> . 
This is a rather modest specific angular momentum, compared 
with observational estimates for the gas in star-forming clouds 
(e.g. Bodenheimer 1995), so the expectation must be that most 
prestellar condensations are flattened by rotation during their 
Second Collapse. 

7.3. Forming close, low-mass binaries by Secondary 
Fragmentation 

If, during The Second Collapse, the spin angular momentum of 
a prestellar condensation is sufficient to cause it to flatten, it is 
likely to fragment into two or more pieces (Tsuribe & Inutsuka, 
1999a, 1999b), and hence to form a binary system, or an un- 
stable higher-order multiple system. We call this Secondary 
Fragmentation. The mean semi-major axis of a binary system 
formed in this way is of order 



G m Wp, 

3 Ty. 



50 AU 



(78) 



5AU 



(79) 



Higher-order multiples are likely to evolve dynamically by 
ejecting lower-mass components, and this may ultimately lead 
to somewhat closer binary systems than predicted by Eqn. (I78> . 

In discussing binary systems with brown-dwarf primaries, 
it is standard practice to include also systems with very low- 
mass hydrogen burning primaries m, ^ 0.1 Mg. These systems 
will therefore have semi-major axes 

(m^ + m,) 
_ 0.1 M, 

Here we have assumed that most of the mass of the orig- 
inal prestellar condensation (the one formed by Primary 
Fragmentation) goes into the two components of the binary, 
with masses m, and . This is a reasonable assumption, given 
that the prestellar condensation is by this stage very tightly 
bound. The inequality in Eqn. il9i acknowledges the possibil- 
ity that the binary is hardened, by dynamical ejection of addi- 
tional components in an unstable higher-order multiple, and/or 
by dissipative interactions with residual gas. 
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Comparing the locus described by Eqn. J78> with the obser- 
vations reported in Close et al. (2003; their Fig. 15), and not- 
ing the good agreement, we speculate that the genesis of close, 
low-mass binaries is due to The Second Collapse. Confirmation 
of this suggestion will require detailed numerical simulations, 
with appropriate initial and boundary conditions, realistic ther- 
modynamics, and physical transport of angular momentum. 
Bate (1998) reports a simulation of Secondary Collapse, which 
fails to result in Secondary Fragmentation. However, this sim- 
ulation uses a barotropic equation of state, rather than a proper 
treatment of the energy equation and the associated radiation 
transport. It may also be influenced by numerical transport of 
angular momentum. We - and others - are currently revisiting 
this problem using an SPH code with radiation transport, but 
conclusive results are not yet available. 

7.4. Forming close, low-mass binaries by Secondary 
Fragmentation in circumstellar discs 

Secondary Fragmentation requires that the primary fragment 
be rotating, and the only Primary Fragmentation scenario in 
which rotation is explicitly considered is disc fragmentation. 
Therefore we now consider whether close low-mass binaries 
can be formed in the outer parts of massive extended discs, and 
whether close low-mass binaries formed in this location may 
subsequently end up in the field. 

We note that a close, low-mass binary formed in this way 
has binding energy 



- Q. 



2a 



^ 6x lO'^^erg 



(m, + m^) 



0.1 M„ 



(80) 
(81) 



where we have assumed equal-mass components (m, = 
and used Eqn. (I79t . For comparison, the binding energy of the 
close, low-mass binary as a whole to the Sun-like star at the 
centre of the original disc is 



G(m, -I- mJM* 
R 

20 X 10^^ erg 



(m, + m^) 



0.1 M„ 



R 



lOOAU 



(82) 



(83) 



It may therefore be possible for tidal encounters occasionally 
to detach a close, low-mass binary from a central Sun-like star, 
without destroying it (i.e. to unbind the wide system without 
unbinding the close one), and thereby to populate the field with 
close, low-mass binaries that were originally formed in circum- 
stellar discs. 

To evaluate this possibility more accurately, we consider a 
perturbing star of mass M^^^^, travelling at speed Vp^^^ and pass- 
ing a binary system with total mass mg,,^ and separation i^,^, 
with the distance of closest approach being D^^,^^. The tide of 
the passing star will deliver a velocity impulse to one compo- 
nent of the binary, relative to the other, given by the product of 
the tidal acceleration and the duration of the interaction. 



Av 



G -^pgRT S^^^ ^PERT 
PERT PERT 



(84) 



strictly speaking, the duration of the interaction (the second 
term on the righthand side of Ean.l84> should be 



MIN 



i/2\ 



but including this factor would only strengthen our conclusion, 
so we omit it in the interests of simplicity. The velocity impulse 
given by Eqn. J84> will unbind the binary if it exceeds 



1/2 



Therefore disruption requires 



m, 



BIN 

.3 



GMl 



D4 v2 

pert pert 



(85) 



(86) 



Using the parameters derived above, we estimate m^^^/s^^^ 
0.5 X 10"^ g cm""* for the typical close, low-mass binary formed 
as a result of The Second Collapse, whereas m^m/^lm ~ ^-^ ^ 
10^'^gcm"-^ for a typical wide binary formed as a result of 
disc fragmentation. Therefore there appears to be a signficiant 
range of perturber parameters which allow the close, low-mass 
system to survive whilst destroying the wide system. 

Thus it is possible that close, low-mass binaries are formed 
in the outer parts of massive circumstellar discs round Sun-like 
stars, in a two-stage process: (i) the formation of a low-mass 
prestellar condensation by Primary Fragmentation (facilitated 
by the fact that the condensation can radiate sufficiently fast 
through thermal dust emission to keep cool); (ii) Secondary 
Fragmentation of this primary fragment (during The Second 
Collapse, i.e. whilst is dissociating) to produce a close, low- 
mass binary system. Additionally, it may sometimes be possi- 
ble to detach a close, low-mass binary formed in this way from 
the Sun-like star at the centre of the natal disc, without destroy- 
ing the close, low-mass binary system - and thereby to populate 
the field with close, low-mass binaries. 

This scenario is attractive because it explains why 
Burgasser et al. (2005) infer that brown dwarfs in wide orbits 
around Sun-like stars have a higher probability of being in a 
close binary system with another brown dwarf than do brown 
dwarfs in the field. However, we must caution that this infer- 
ence is presently based on small-number statistics and needs to 
be confirmed with a larger sample. 

8. Conclusions 

We have reviewed the thermodynamic processes which are pre- 
sumed to determine the minimum mass for star formation - 
namely (i) Primary Fragmentation and the Opacity Limit, (ii) 
Secondary Fragmentation during H, dissociation - and the con- 
sequences these thermodynamic processes have for the statis- 
tics of binary systems containing brown dwarfs. 

8.1. Primary Fragmentation and the Opacity Limit 

Here the presumption is that a prestellar condensation only 
fragments if it can radiate fast enough to stay approximately 



A. P. Whitworth and D. Stamatellos: The minimum mass for star formation, and the origin of binary brown dwarfs 



13 



isothermal. We have treated three different generic star forma- 
tion scenarios, (a) hierarchical fragmentation of a 3-D medium, 
(b) one-shot 2-D fragmentation of a shock-compressed layer, 
and (c) fragmentation of a circumstellar disc. We believe that 
these three scenarios cover the different basic situations in 
which minimum-mass fragments may be created (as regards 
geometry, background dynamics, competing processes, etc.). 

Hierarchical 3-D fragmentation represents fragmentation in 
a cloud undergoing overall collapse, with no preferred direc- 
tion (i.e. a statistically isotropic situation). The main processes 
competing (successfully) against fragmentation are contraction 
of the background and accretion. We re-derive the expression 
for the minimum mass obtained by Rees (1976) under the as- 
sumption that the fragment is marginally optically thick (Eqn. 
\13i . and demonstrate that this assumption is valid in the lo- 
cal contemporary interstellar medium. We also derive a simple 
analytic expression for the luminosity of a spherical fragment 
which is not necessarily marginally optically thick (Eqns.|9]& 
not , and hence we obtain expressions for the minimum mass 
which can be used in circumstances other than contemporary 
local star formation, when it cannot necessarily be assumed a 
priori that the minimum-mass fragment is marginally optically 
thick (Eans.[T4l& ll5t . We also rehearse the reasons why hier- 
archical 3-D fragmentation may not be a useful paradigm for 
contemporary star formation. 

One-shot 2-D fragmentation represents the situation where 
star formation is triggered by two (or more) supersonically col- 
liding streams, producing a shock-compressed layer (or fila- 
ment). Therefore (a) the geometry is anisotropic, (b) ram pres- 
sure is important, and (c) the fragments are initially flattened. 
Here the main processes competing against fragmentation - 
and thereby fixing the scale of fragmentation - are the con- 
tinuing inflow of material into the layer, and the extra heating 
associated with the resulting accretion shock. The masses of 
fragments depend on the flux of matter in the colliding flows 
(Eqn. I39> . and the initial aspect ratio of a fragment is approxi- 
mately equal to the Mach Number of the accretion shock (Eqn. 
I33I I. Expressions are derived for the time at which the accu- 
mulating layer fragments (Eqn. I29> . for the properties of the 
resulting fragments (Eans. l30lto I32>. for the luminosity of a 
flattened disc-shaped fragment (Eans.l37l&l38>. and hence for 
the minimum mass in situations where the fragment is opti- 
cally thin, marginally optically thick, or strongly optically thick 
('Eans.l?Ulto l42l i. We also formulate the condition which must 
be satisfied for minimum-mass fragments to be optically thin or 
optically thick. Hence we show that minimum-mass fragments 
in the outer galaxy should be optically thin at their inception, 
and conversely minimum-mass fragments in the inner galaxy 
should be optically thick at their inception. 

Disc fragmentation represents the situation where the star 
forming material is in orbit about an existing star and therefore 
(a) the geometry is again anisotropic, (b) centrifugal accelera- 
tion is important. The main process competing with condensa- 
tion of a proto-fragment is tidal shear, which tends to tear proto- 
fragments apart, unless they can cool and condense out on a 
dynamical timescale. We present a derivation of the Toomre 
criterion which has the advantage both of being very simple 
(Eans.l47l&l48t. and of being appHcable to non-equihbrium 



discs. This is important because in real star formation events 
discs may not have sufficient time to settle into equilibrium be- 
fore they fragment. In a disc, the masses of fragments depend 
on the local sound speed, a, and the local vorticity, w, since 
condensation of a fragment is resisted by both pressure and ro- 
tation. Hence the minimum mass depends on the properties of 
the central star (mass, M*, and luminosity, L*), but with rea- 
sonable assumptions this dependence is quite weak. 

The fact that, for contemporary local star formation, all 
three generic scenarios predict a minimum mass in the narrow 
range 0.001 to 0.004 M^, is remarkable. In the different scenar- 
ios, the expressions for the Jeans mass are completely differ- 
ent, the geometries are different, the competing processes are 
different, and the radiation transport regimes are different. In 
hierarchical 3-D fragmentation, OTjeansi depends only on p and 
a, and the shortest condensation time is for the largest frag- 
ment. In one-shot 2-D fragmentation of a shock-compressed 
layer, /Wje^nsz depends also on v and f, fragments are initially 
flattened, the critical quantity is the flux of matter flowing into 
the layer, and the shortest condensation time is for a finite-mass 
fragment. In a disc, ^Jjeansd depends on a and Z, hence on M*, 
Li, and R, the critical quantity is aP'/a^, and the minimum-mass 
fragment is optically thick (f ~ 100), unlike in the other two 
scenarios (where f ~ 1). Yet despite these differences, the min- 
imum mass has a very small range. This implies that in contem- 
porary, local star formation, the details of geometry and back- 
ground dynamics have little influence on the minimum mass. 
(It might be tempting, but we believe it would be incorrect, to 
assert that this result does not need proof.) 

In addition to evaluating the minimum mass, we have 
shown that it is hard for discs to undergo gravitational fragmen- 
tation at small radii, basically because the timescale on which a 
proto-fragment is sheared apart is shorter than the timescale on 
which it can cool and condense out. Consequently there should 
be a dearth of brown dwarfs in orbit around Sun-like stars, 
and this provides a possible explanation for the observationally 
inferred Brown Dwarf Desert. Conversely, brown dwarfs can 
condense out at large radii in massive extended discs, because 
here the timescale for a proto-fragment to be sheared apart is 
likely to exceed the cooling timescale. We therefore suggest 
that disc fragmentation is the most likely formation mechanism 
for those brown dwarfs seen in wide orbits about Sun-like stars. 
We also point out that fragmentation of the outer parts of mas- 
sive extended discs should occur very quickly, essentially on a 
dynamical timescale (~ lO'* years), and therefore the paucity of 
massive extended discs need not mean that they do not form, 
but simply that they quickly self-destruct by condensing into 
brown dwarfs and low-mass H-burning stars. 

8.2. Secondary Fragmentation during dissociation 

A prestellar condensation formed by Primary Fragmentation 
can fragment further during the Second Collapse, when disso- 
ciation of H, acts as a sink for the internal energy being de- 
livered by compression. If at this stage the protostellar con- 
densation is flattened by rotation, the Secondary Collapse can 
lead to Secondaiy Fragmentation and the formation of a close 
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low-mass binary system. From simple thermodynamic argu- 
ments, H, dissociation occurs at T^^^^ ~ 2000 K and p^^^^ ~ 
10"' g cm Therefore the resulting binary separations should 
approximate to 



fl ~ 5AU 



0.1 Mo 



(87) 



If this is compared with the plot of Close et al. (2003; their Fig. 
15), it is seen to provide a remarkably good fit to the observed 
separations of close low-mass systems, suggesting strongly that 
these systems have been formed by Secondary Fragmentation. 

Since Secondary Fragmentation requires that the primary 
fragment be rotating, and since the only fragmentation scenario 
we have considered which explicitly takes account of the spin 
of a proto-fragment is disc fragmentation, it is tempting to con- 
sider the outer parts of discs as the location for forming close 
low-mass binary systems. This suggestion derives support from 
the observation that brown dwarfs in wide orbits about Sun-like 
stars appear more likely to have brown-dwarf companions than 
brown dwarfs in the field (Burgasser et al. 2005). The simplest 
explanation for this would seem to be that low-mass prestel- 
lar condensations form by Primary Fragmentation in the outer 
parts of discs (as described in Section |6} and then undergo 
Secondary Fragmentation to form close low-mass binaries, in 
situ. 

Furthermore we show that close low-mass binaries formed 
in this way can sometimes survive being impulsively separated 
from the Sun-like star, and therefore brown-dwarf binaries in 
the field could also have formed by disc fragmentation. 

Appendix A: Radiation transport in a one-zone 
fragment 

A.1. A spherical fragment 

To justify both Eqns. and ilOi . we consider a uniform- 
density spherical fragment of mass m and radius r, and treat 
separately the optically thin and optically thick limits. 

In the optically thin limit, the term f^(r) can be neglected 
in comparison with f '(T), and so Eqn. (|9j reduces to 



Ahin = 4mK^{T)o-,^ T 



(A.l) 



This is an exact expression, in the sense that it defines the 
Planck-mean opacity, Kp{T). 

In the optically thick limit, the term f^^\T) can be ne- 
glected in comparison with f^(T), and so Eqn. now reduces 
to 



r 



16 71 T'^ 



(A.2) 



Since f^{T) is the Rosseland-mean optical-depth between the 
centre of the spherical fragment and its surface, we can define 
the diffusion time. 



(A3) 



which is just the time it takes radiation to random-walk from 
the interior of the spherical fragment to its exterior. Ehminating 



T^(T) in favour of fp,pp, and putting cr^^ - a^^c/A- (where a^^ is 
the radiation energy-density constant), Eqn. (IA.2> becomes 



r 

-*-'TI 



(A.4) 



The righthand side of Eqn. (IA.4> is the total radiant energy in- 
side the spherical fragment divided by the time this radiation 
takes to escape. In the optically thick limit this is a reasonable 
approximation to the luminosity. 

Finally we note that, when the spherical fragment is neither 
optically thin, nor optically thick, and fj,(r) ^ f^{T) ^ 1, the 
luminosity of the fragment is given by 

2 

3 ' 

so it cools almost as well as a blackbody. 



(A.5) 



A.2. A disc-like fragment 

To justify both Eqns. ( I37> and ( I38> — and Eqns. i61\ and (|62 



with r, — » r^jscD — we consider a uniform-density disc-like 
fragment of mass m,, radius r,, and thickness 2zj (where the 
assumption is that z r, i.e. the disc is geometrically thin). 
We again treat separately the optically thin and optically thick 
limits. 

In the optically thin limit, the term fj^(T) can be neglected 
in comparison with fp"'(r), and so Eqn. j37> reduces to the 
exact form: 

= 4m^_-K^(T)cr^^T\ (A.6) 

In the optically thick limit, the term '(T) can be ne- 
glected in comparison with fj^{T), and so Eqn. i31l now re- 
duces to 

4;rr2(r,^ 

The time for radiation to diffuse from the interior of the disc- 
like fragment to its flat surfaces at +z is 

W ^ (A.8) 

where the factor 2 arises because the radiation must diffuse 
in the z-direction to escape. It follows that Eqn. iA.7\ can be 
rewritten as 

-^THICK ~ T ■ (A. 9) 

^DIFF 

The righthand side of Eqn. ( IA.9> is the total radiant energy in- 
side the disc-like fragment, divided by the time this radiation 
takes to escape, which in the optically thick limit is a reason- 
able approximation to the luminosity. 

When the disc-like fragment is neither optically thin, nor 
optically thick, and f^iT) ^ fp{T) ^ 1, its luminosity is given 
by 

X^^, ^ 27rrV,, X 1 , (A. 10) 

so it cools exactly like a black body from its flat surfaces. 
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